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ABSTRACT 

ASCA X-ray spectra of many ULXs were described in terms of optically thick emission from hot {kTm ~ 1 — 2 keV) 
accretion disks, while recent XMM-Newton and Chandra observations have revealed a cool {kTin ^ 0.2 keV), soft X-ray 
excess emission from a number of them. Here we utilize improved calibration and high signal-to-noise XMM-Newton 
spectra of NGC 1313 X-1 to present evidence for a cool (~ 0.2 keV) soft excess and a curved or a cutoff power-law 
(F ~ 1 — 1.5, Ecutoff 3 — 8 keV). The high energy curvature may also be described by a hot 1 — 2.5 keV) 
multicolor disk blackbody. The soft excess emission is unlikely to arise from a cool disk as its blackbody temperature 
is similar in three XMM-Newton observations, despite a change in the observed flux by a factor of about two. Thus, 
previous estimates of the black hole mass of lOOOM© for NGC 1313 X-1 based on the temperature of the soft excess 
emission is unlikely to be correct. The power-law cutoff energy is found to decrease from ~ 8 keV to 3 keV when the 
ULX brightened by a factor of about two. The unusual spectral properties of NGC 1313 X-1 are difficult to understand 
in the framework of the disk/corona models generally adopted for the black hole binaries or active galactic nuclei and 
may require to invoke super-critical accretion rates. 

Subject headings: accretion, accretion disks — stars: individual (NGC 1313 X-1) — X-rays: stars 



1. INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are point X-ray 
sources with luminosities exceeding ~ 10'^^ erg s~^, out- 
side the nucleus of a galaxy. The nature of ULXs contin- 
ues to be an enigma, since their adopted isotropic high 
energy output surpasses the Eddington limit of even the 
most massive stellar mass black holes (BHs), sometimes 
by large factors. Several models have been proposed to 
explain the high luminosities of ULXs. The most popu- 
lar is the "intermediate mass black hole (IMBH)" with 
mass Mbh ^ 10^ - 10"* M© (e.g., Colbert & Mushotzky 
1999, hereafter CM99). Other models include (i) XRBs 
with anisotropic emission (King et al. 2001), (ii) beamed 
XRBs with rclativistic jets directly pointing towards us i. 
e., scaled down versions of blazars (Mirabel & Rodriguez 
1999), and (iii) XRBs with super-Eddington accretion 
rates (Begelman 2002). 

Several observations suggest that ULXs may be simi- 
lar or scaled up versions of the Galactic X-ray binaries 
(XRBs). Orbital modulations in some ULXs (Bauer et 
al. 2001; Sugiho et al. 2001) imply their binary na- 
ture. Observations of spectral transitions between the 
low (hard) and the high (soft) state from two ULXs in 
IC 342 (Kubota et al. 2001) demonstrate their similarity 
to the Galactic BH binaries. The ASCA X-ray spectra 
of ULXs have been described as the emission from opti- 
cally thick accretion disks (Makishima et al. 2000 and 
references therein). However, the inferred inner disk tem- 
peratures in the range of ~ 1 — 2 keV are too high for an 
ULX accreting at sub-Eddington rates. Recent Chandra 
and XMM-Newton observations of ULXs have shown soft 
X-ray excess emission which has been interpreted as the 
optically thick emission from thin accretion disks with 
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temperatures in the range of ^ 100 — 300 eV, suggesting 
intermediate mass black holes (~ 10^ - lO'' Mq) accret- 
ing at sub-Eddington (~ OALEdd) rates (see e.g.. Miller, 
Fabian & Miller 2004a,b; Miller & Colbert 2003). Thus 
ULXs may either be ~ 20 — IOOMq black holes accreting 
at near or super Eddington rates or ~ 1000 — 5OOOM0 
black holes with sub Eddington rates. 

One of the first ULXs observed by XMM-Newton to 
show the soft excess emission that led to the ~ lOOOM© 
IMBH interpretation, is NGC 1313 X-1 (Miller et al. 
2003). Here we revisit the spectral analysis of the ULX 
NGC 1313 X-1 using aU pubhcly available XMM-Newton 
observations with improved calibration. We show that 
the high energy continuum of the ULX is curved and can 
be well described either by a cutoff power law or a high 
temperature accretion disk component. We also confirm 
the soft excess emission earlier detected by Miller et al. 
(2003). 

2. OBSERVATION AND DATA REDUCTION 

The ULX NGC 1313 X-1 has been observed four- 
teen times by XMM-Newton with exposure times rang- 
ing from 9 ks to 41 ks. We used the EPIC data obtained 
from the nine observations that are publicly available and 
listed in Tabled The first XMM-Newton observation on 
2000 October 17 was analyzed by Miller et al. (2003) 
and led to one of the first discoveries of a cool soft excess 
component from an ULX, and a BH mass estimate of 
~ IO^Mq. They used the MOS data alone, as the PN 
data could not be processed with the SAS 5.3.3 avail- 
able in 2002. There were also significant differences in 
the PN and MOS spectra below - 0.6 keV: the MOS 
data showed generally flatter spectra compared with the 
PN data. Here we used SAS 6.5.0 with updated cal- 
ibration to process and ffiter both PN and MOS data 
obtained from all the observations listed in Tabled All 
the observations were affected by a high particle back- 
ground. Cleaning of the flaring particle background re- 
sulted in reduced 'good' exposure times listed in Tabled 
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TABLE 1 

Log of publicly available XMM-Newton observations of 
NGC 1313. 



Observation 


Date of 


Usable 


ID 


Observation 


exposure 


0106860101 


2000-10-17 


23.5 ks 


0150280101 


2003-11-25 


2.5 ks 


0150280201 


2003-12-09 


0.0 ks 


0150280301 


2003-12-21 


7.4 ks 


0150280401 


2003-12-23 


3.2 ks 


0150280501 


2003-12-25 


1.7 ks 


0150280601 


2004-01-08 


6.2 ks 


0150280701 


2003-12-27 


0.0 ks 


0150281101 


2004-01-16 


3.5 ks 



Events in the bad pixels and those in adjacent pixels 
were discarded. Only events with pattern — 4 (single 
and double) for the PN and 0-12 for the MOS were 
selected. Only three observations, with listed exposure 
times > 6 ks, are useful for our purpose of detailed spec- 
tral analysis. 

3. ANALYSIS & RESULTS 

We extracted PN and MOS spectra from all the obser- 
vations using circular regions with radii of 35" centered 
at the position of NGC 1313 X-1. We also extracted PN 
and MOS background spectra using appropriate nearby 
circular regions free of sources. We created spectral re- 
sponse files using the SAS tasks rmfgen and arfgen. 
Spectral bins were chosen such that there was a mini- 
mum of 20 and 15 counts per spectral channel for the 
PN data of the first observation and all other data sets, 
respectively. These spectra were analyzed with XSPEC 
11.3. The errors on the best-fit spectral parameters are 
quoted at a 90% confidence level. 

First we fitted a simple absorbed power law (PL) model 
to the PN and MOS spectra of NGC 1313 X-1 obtained 
from the first observation. We fit the MOSl and M0S2 
data jointly with an overall normalization constant to ac- 
count for possible differences in source extraction areas 
or calibration uncertainties. We used the 0.2 — 12 keV 
band in all the fits. The simple power law model re- 
sulted in minimum — 749.3 for 613 degrees of free- 
dom (dof) and 712.1 for 656 dof for the PN and MOS 
data respectively, thus providing statistically unaccept- 
able fits to both the PN and MOS data. We have plot- 
ted the ratios of the PN data and the best-fit PL model 
in Figure and the ratios of the MOS data and the 
best-fit model in Fig. ^p. Both the plots are similar and 
clearly show a broad hump in the 3 — 8 keV band and 
a relatively narrow soft hump or excess emission below 
1 keV. Addition of a soft (~ 0.1 keV) blackbody (BB) 
component improves the fit to the PN data significantly 
(Ax^ = —137.8 for two additional parameters). Similar 
improvements are found for the MOS data also. Fig. QJ; 
shows the ratio of the PN data and the best-fit BB+PL 
model. The broad hump is now seen as excess emission 
in the ~ 2 — 7 keV band and in the lack of emission 
above ~ 7 keV. This change in the shape of the broad 
hump from Fig. ^ to Fig. ^ is due to the flattening 
of the power-law with the addition of the soft black- 
body component. Addition of a multicolor accretion disk 
blackbody (MCD; diskbb in XSPEC) further improves 



the fit to the PN data (Ax^ = 16.7 for two additional 
parameters). This is an improvement at a significance 
level of > 99.99% based on the maximum likelihood ratio 
(MLR) test. To further verify the high energy turnover 
of the X-ray spectra, we replaced the PL component in 
the BB-I-PL model by a cut-off PL. The BB plus cut-off 
PL model yielded Ax^ = —17.10 for one additional pa- 
rameter as compared to the BB-I-PL model. This is an 
improvement at a significance level of > 99.99% based 
on the MLR test. Thus we conclude that the presence of 
the broad hump emission is statistically required. 

Since we did not notice significant differences in the PN 
and MOS data, we present results based on joint spectral 
fitting of the PN and MOS data. TablcHlists the best-fit 
parameters for the spectral models BB+PL, BB+cutoff 
PL, BB-^MCD and BB+MCD+PL. Good signal-to-noise 
spectra are required in order to infer the presence of the 
broad hump or high energy turnover. We used only three 
observations in which the ULX was detected at a level 
of 3ct or better in the 8 — 10 keV band. The X-ray spec- 
tra of NGC 1313 X-1 obtained on 2000 October 17 and 
2003 December 21 clearly statistically require a cutoff in 
the power law or a hot MCD component in addition to 
the cool soft component, while the lower signal-to-noise 
spectra of 2004 January 8 are consistent but do not sta- 
tistically require the presence of the cutoff or the MCD 
component. 

We show the spectral variability of NGC 1313 X-1 in 
FigureElin terms of the ratios of the data and the best-fit 
simple absorbed power-law model derived from a joint fit 
to the PN data obtained from the observations of 2000- 
10-17 and 2003-12-21. In the joint fit, the power-law in- 
dex and the normalization were allowed to vary indepen- 
dently while the absorption columns for the two spectra 
were tied to each other and varied jointly. The best-fit 
column is ~ 2.4 x 10^^ cm~^ and the photon indices are 
similar (F ~ 2.1) for both the spectra. The ULX was 
about a factor of two brighter in the 2003-12-21 obser- 
vation. There is a clear decrease in the power-law cutoff 
energy (see Fig. Eland Table EJ. 

4. DISCUSSION 

Based on the spectral analysis of the available high 
signal-to-noise PN and MOS data and using the im- 
proved calibration of the EPIC instruments, we found 
that the X-ray spectrum of NGC 1313 X-1 is curved at 
high energies > 6 keV and shows a soft excess emission 
below 2 keV. This ULX showed spectral softening with 
increasing flux, similar to that observed from BH XRBs 
and Seyfert 1 galaxies. The 0.2 — 10 keV luminosity of 
the ULX was 4.3 x lO^^ erg s^^ in 2000 October 17 with 
an increase to 7.9 x 10'^^ erg s~^ in 2003 December 21, 
assuming a distance of 3.7 Mpc (TuUy 1988). 

4.1. A cool corona 

The soft excess and high energy curvature in the X-ray 
spectra of NGC 1313 X-1 are well described by a model 
comprising a cool BB (kT ~ 0.2 keV) and a PL (F ^ 
1 — 1.5) with a high energy cutoff, Ecutoff ~ 3 — 8 keV 
(see Table|21l. The soft excess and power- law, with a high 
energy cutoff, are characteristic spectral components ex- 
pected from the accretion disk-corona systems in BH 
XRBs and AGNs. The soft excess is thought to be the op- 
tically thick emission from an accretion disk. The cutoff 
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Fig. 1. — The ratios of observed data obtained on 2000 October and the best-fit models. The data are binned heavily for display purposes 
only. Ratios of (a) observed PN data and the best-fit absorbed PL model, (b) observed MOS data and the best-fit absorbed PL model, 
(c) the PN data and the best-fit absorbed BB-I-PL model, and (d) the PN data and the best-fit absorbed BB-I-MCD+PL model. A clear 
cutoff at 7 keV is seen in panel (c). 
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Fig. 2. — Ratios of the data and the best-fit absorbed power-law 
model jointly fitting the two PN spectra of 2000-10-17 and 2003- 
12-21. The absorption columns for the two spectra were tied to 
each other, while the power law indices and normalizations were 
independently. The photon indices are similar (F ~ 2.1) for both 
the spectra, despite a factor of two increase in the normalization in 
the observation of 2003 December 21. It is clear that the power-law 
cuts off at lower energy in the 2003 observation. 



power-law component arises from thermal Comptoniza- 



tion with a cutoff energy characteristic of the tempera- 
ture of the hot comptonizing corona. The observed pho- 
ton index of NGC 1313 X-1 is similar to that of BH XRBs 
in a low/hard state. In many such BH XRBs and also in 
Seyfert galaxies, the cutoff energy is usually observed in 
the range of ^ 100 keV (see e.g., McClintock & Remil- 
lard (2004). If the soft excess and the cutoff power-law 
emission from NGC 1313 X-1 arise from a disk-corona, 
then the corona must be extremely cool compared to the 
coronae of both XRB and AGNs. This is a rather puz- 
zling result. There are also problems with the interpre- 
tation of the soft excess emission as the optically thick 
emission from a thin disk (see below). 

4.2. The soft X-ray excess emission 

The temperature of the cool soft component, parame- 
terized as a blackbody, is ^ 0.2 keV in all three observa- 
tions. The blackbody normalization increased by a factor 
of ~ 2 - 3 from 2000 October 17 to 2003 December 21. 
The cool soft emission was the weakest amongst the three 
spectral components in the observation of 2000 October 
17, contributing only 17% to the 0.2 - 10 keV emis- 
sion. In the observations of 2003 December 21 and 2004 
January 8, the cool soft component contributed ~ 30% 
and 6%, respectively, in the 0.2 — 10 keV band. Miller 
et al. (2003) interpreted this component as the optically 
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thick, thermal emission from a cool accretion disk and 
estimated a black hole mass of ~ 10"^ Af© for NGC 1313 
X-1. 

Roberts et al. (2005) have pointed out a number of 
problems with the above interpretation, which is based 
on the extrapolation of the X-ray spectra of XRBs to 
the IMBH regime and the detection of soft X-ray ex- 
cess emission from ULXs. It is also assumed that the 
accretion disks of ULXs extend to the last stable or- 
bit, i.e., they are in a 'high' state. This corresponds 
to the high/soft state of BH XRBs, where the optically 
thick, thermal emission from the disk completely domi- 
nates the 2 — 10 keV spectrum. This is generally not the 
case in ULXs, where the cool soft X-ray excess emission 
only contributes 10 — 30% to the total X-ray emis- 
sion. The MCD component in BH XRB is associated 
with a steep power law (photon index F > 2.4); this is 
not always the case in ULXs (see Miller et al. 2003, 
2004; Roberts et al. 2005; Dewangan et al. 2005). The 
BB+PL model resulted in a photon index of 1.74 ± 0.04 
for NGC 1313 X-1 in the first XMM-Newton observation 
(see Table 111 . Furthermore, the disk temperature of BH 
XRBs increases with their luminosity (see e.g.. Miller, 
Fabian & Miller 2004), but the temperature of the soft 
component from NGC 1313 X-1 is remarkably constant 
in spite of the ~ 80% change in the 0.2 — 10 keV lumi- 
nosity. There is also a definite curvature or cutoff above 
- 6 keV (2000 October) or ~ 3 keV (2003 December) 
in the X-ray spectrum of NGC 1313 X-1. This curva- 
ture is usually not present in the high state of XRBs: 
instead, a steep power-law (F 2.5 — 4) extending up to 
several hundred keV is observed. GRS 1915-1-105 and 
GRO J1655-40 are the XRBs that show e-folding cutoffs 
at ~ 3.5 keV in their high state (McClintock & Remil- 
lard 2004), but, unhke NGC 1313 X-1, their power-law 
slopes are very steep. 

What is the origin of the soft excess component if it 
is not the thermal emission from a cool accretion disk? 
There are other classes of accreting objects that show 
similar soft components. Many accreting X-ray pul- 
sars also show soft X-ray excess emission with a black- 
body temperature of 0.1 - 0.2 keV (Paul et al. 2002; 
Hickox et al. 2004; Neilsen et al. 2004). Narrow-line 
Seyfert 1 galaxies (NLSls) and quasars show soft excess 
emission below 1 keV that is well described by a black- 
body. The apparent temperature of the soft component, 
kT ~ 0.1 — 0.3 keV, is similar in a diverse selection of 
active galactic nuclei (e.g., Czerny et al. 2003; Gierlihski 
& Done 2004). This observation can be compared with 
the similarity of the temperature of the soft component 
of NGC 1313 X-1 at different luminosities. The ori- 
gin of the soft excess emission is not well understood. 
King & Pounds (2003) suggest that BHs, accreting at or 
above the Eddington rate, produce winds that are opti- 
cally thick in continuum. Blackbody emission from these 
winds may provide the observed soft excess emission from 
ULXs. 

4.3. Comparison with GRS 1915+105 

XMM-Newton observed the Galactic micro-quasar 
GRS 1915-1-105 in its "radio-loud plateau state" (Martoc- 
chia et al. 2005). There is a very striking similarity be- 
tween GRS 1915+105 in its plateau state and NGC 1313 
X-1 observed on 2000 October. Both show small ampli- 



TABLE 2 

Best-fit spectral model parameters for NGC 1313 X-1. 



Model Date of Observation 

Parameters 2000-10-17 2003-12-21* 2004-01-08 

BB+PL 





2 44+"- *4 


o 7+0. b 
' -0.4 


3 0+0'^ 


kTggi keV) 


16+""^ 


14+0-04 


n -,0+0.04 


Abb"" 


1 1+0.3 
^•^-0.2 


2 1+4-5 
^-1-1.5 


1 0+0-4 

1-^-0.6 




.0.85 


1.6 


0.85 


r 


1 74+0.04 


9 40+0.08 


2 1+0-1 




4-0-0.3 




12.4^2-2 




3.45 


13.0 


7.5 


fx'' 


2.6 


4.8 


4.1 


Lx" 


4.3 


7.9 


6.7 




1250.6/1271 


353.0/372 


534.6/608 


BB-l-cutofr PL 


AfffClO^i cm-2) 


^-^-0.3 


9 O+0.3 


9 7+0.5 
^-'-0.3 


kTBBi keV) 




19+0-04 


20+0-04 


Abb"" 


1 0+''-^ 


9 O+0.7 

^-■'-o.s 


-1 9+0.5 


Ibb"" 


1.0 


1.8 


1.0 


r 


1 91+0.12 
^-^^-0.29 


i.38i?:l^ 


l-<5O-0.40 


Ecutoffi keV) 




3.6l}-2 


> 0.36 


Apl" 


4Q+0.2 




11 3+2-1 
11-^-2.6 




2.8 


6.8 


6.3 


fx" 


2.6 


4.8 


4.0 


Lx" 


4.3 


7.9 


6.5 


X™,,;„/dof 


1231.7/1270 


337.2/371 


533.6/645 


BB+MCD 


AfffClO^i cm-2) 




2 0+0-y 


1 5+0-2 
1-3-0.2 


fcTss( keV) 




24+0-02 
'^•^^-0.02 


0.28+ 


Abb"" 


1 Q + O.Ol 


9 7+0.6 
^- '0.4 


1 9+0.2 
1-^-0.2 


Ibb"" 


0.85 


2.2 


1.6 


kTMCoi keV) 


r, r,t- + 0.09 

^-^O_0.11 


1 Qg-fO.14 


-, OO+0.20 
l-»»-0.17 


Amcd'' 


4 7+1-0 
^- '-0.6 


is.otti 


Imcd" 


2.3 








2.6 


4.8 


3.3 


Lx" 


4.3 


7.9 


5.4 


xi,„/dof 


1259.6/1271 


344.3/372 


547.1/608 


BB-I-MCD+PL 


Afj/(102i cm-2) 


01-1 + 0.53 
Z.OU_Q 32 


9 4+0.5 


9 O+0.9 


kTBB{ keV) 


n 10+0.02 


21+0-03 


n iq+o.o5 


Abb"" 


1 o+°-2 

^■"-0.2 


9 7+0.8 
^- '-0.5 


1-1-0.7 


Ibb'' 


0.80 


2.2 


0.9 


kTMCoi keV) 


9 cc+O.SO 


1-1-0.2 


1-^-1.9 


Amcd'' 


1-3-0.8 


1S4 -1+104.3 
l'54-l-64.6 


9 n + 11.1 
■^■"-2.0 




1.2 


3.8 


0.5 


r 




1 4+1-0 


9 9+0.8 
^■^-1.0 


Apl" 


3.511- 






f b 
JPL 


2.2 


1.4 


6.8 




2.6 


4.9 


4.0 


Lx" 


4.3 


8.0 


6.5 


xi,./dof 


1231.5/1269 


335.8/370 


534.0/606 



''BB normalization in units of 10 ^ X -pjj^p-i where L39 is the 

source luminosity in units 103^ erg s~*^ and DIO is the distance to 
source in units of 10 kpc. 
iTntrinsic flux in units of 10^*^2 ergs cm~2 s~i and in the 0.2 — 
10 keV band. 

'^MCD normalization in units of 10~3( '^'^^g'" )2 x cos9, where Tin 
is the inner disk radius and 6 is the disk inclination angle. 

'l Power-law normalization in units of 

104 photons keV-l cm-2 s"! at 1 keV. 

•^Observed flux and luminosity in units of 10~*^2 ergs cm~2 s~*^ 
and 103^ erg s~^, respectively, and in the band of 0.3 — 10 keV. 

^Source in MOS chip gaps, only PN data were used. 
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tude variability. Simple absorbed power-law fits to the 
X-ray spectra of both the objects result is similar photon 
index (F ~ 1.7), the residuals show similar deficits above 
8 keV and soft excess emissions around 1 keV. Fig. ^ 
is quite similar to Figure 5 in Martocchia et al. (2005). 
The above similarity in the X-ray emission of the XRB 
and the ULX is observed at a factor of ^ 10 difference 
in luminosities. Thus if the X-ray spectral state is solely 
driven by accretion rate relative to the Eddington rate, 
the BH mass of NGC 1313 X-1 is likely a factor of 10 
higher (- 140Mq) than that of GRS 1915+105. The 
deficit and the soft excess emission in the spectrum of 
GRS 1915-1-105 have been interpreted as the reprocess- 
ing from an ionized disk and reflection from optically 
thin, disk wind, respectively. There are minor differ- 
ences which, if real, may lead to the different interpreta- 
tion of the X-ray spectra and hence invalidate the direct 
scaling of their BH masses. The soft excess emission of 
GRS 1915+105 is not well described by a blackbody. In 
addition, GRS 1915+105 shows soft X-ray emission lines 
as well as iron Ka line that are not seen in the spectrum 
of NGC 1313 X-1. The ionized disk reflection cannot ex- 
plain the cutoff or deficit above 3 keV in the 2003-12-21 
spectrum. 

4.4. A hot accretion disk 

The X-ray spectra of NGC 1313 X-1 are also well de- 
scribed by the superposition of a cool BB, a hot MCD 
and a power-law. However, the temperature of the MCD 
component (fcT ~ 2.5 keV on 2000 October) is too high 
for a thin standard accretion disk. This is similar to 
the "high temperature" problem of ULXs observed with 
ASCA (Mizuno et al. 1999; Colbert & Mushotzky 1999; 
Makishima et al. 2000). One possible explanation is in 
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terms of a thin disk around a Kerr BH, since the inner 
radius of a Kerr disk can be = 6 times smaller than that 
of a Swarzschild disk, and the temperature of a thin disk 

varies as r^^^'^ , thus resulting in high inner disk temper- 
atures. Ebisawa et al. (2003) showed that Kerr models 
work well for ULXs but the data imply that the disks are 
highly inclined (> 80°). This model cannot also explain 
the soft excess emission below 2 keV from NGC 1313 
X-1. 

Another attractive scenario is the slim disk model, in 
which the disk luminosity can exceed the Eddington lu- 
minosity by up to a factor of 10 (Abramowicz et al. 
1988, Ebisawa et al. 2003; Kawaguchi 2003). The slim 
disk is geometrically thick and can be much hotter than 
the standard disk. If this is the case, NGC 1313 X-1 must 
be radiating at super-Eddington rates (~ 1 — lOLEdd)- 
Its bolometric luminosity, L^oi ^ 1-5 x 10""^ erg s~^ esti- 
mated using the BB+cutoff PL model for the 2003 De- 
cember observation, implies a BH mass in the range of 
^ 12 — 125M0. A natural consequence of the slim disk 
model is that the soft excess emission can be produced in 
the optically thick winds resulting from the super-critical 
accretion rates as discussed above (Mukai et al. 2003; 
King & Pounds 2003). 
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